Electron Transport in a Multi-Channel One-Dimensional Conductor: Molybdenum 

Selenide Nanowires 
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We have measured electron transport in small bundles of identical conducting Molybdenum Se- 
lenide nanowires where the number of weakly interacting one- dimensional chains ranges from 1-300. 
The linear conductance and current in these nanowires exhibit a power-law dependence on tem- 
perature and bias voltage respectively. The exponents governing these power laws decrease as the 
number of conducting channels increase. These exponents can be related to the electron-electron 
interaction parameter for transport in multi-channel 1-D systems with a few defects. 

PACS numbers: 72.15.Nj,73.23.Hk,73.90.f 
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Interacting electrons in one-dimensional (ID) metals 
constitute a Luttinger Liquid (LL)Q, in contrast to a 
Fermi liquid (FL) in 3-dimcnsional (3D) metals. Trans- 
port properties of ID conductors are strongly modified 
as adding an electron to a ID metal requires changing 
the many-body state of its collective excitations. This 
results in vanishing electron tunneling density of states 
at low energy. Power-law dependent suppression in tun- 
neling conductance has been observed in many systems, 
including fractional quantum Hall edge states |2j, sin- 
gle and multi-walled carbon nanotubes 0, 0, U bun- 
dles of NbSe 3 nanowires [|| and conducting polymers • 
A cross-over from a truly ID Luttinger-liquid (LL) to 
a 3D Fermi-liquid (FL) is expected as ID conductors 
are coupled together, increasing the number of weakly 
interacting channels 0, 0. This transition, however, 
has not been observed in the above (quasi) 1-D systems 
due to the experimental difficulty in preparing identi- 
cal conducting quantum wires to form conductors with 
a few weakly interacting channels. In this letter, we re- 
port temperature and bias dependent electric transport 
measurements on small bundles of Molybdenum Selenide 
(MoSe) nanowires 0, fill H^ |. whose diameter ranges 
from 1-15 nm. These nanowires, which consist of bundles 
of weakly interacting and electrically identical ID MoSe 
molecular chains, show a power-law dependent tunnel- 
ing conductance. The exponent governing the power-law 
decreases as the bundle diameter increases, indicating 
a transition from ID to bulk transport with increasing 
number of conducting channels. 

Crystalline bundles of MoSe chains are obtained from 
the dissolution of quasi-lD Li2MoeSe6 crystals in polar 
solvents. Single crystal Li2Mo6Se6 was prepared as de- 
scribed previously [lj| . X-ray diffraction analysis showed 
hexagonally close packed molecular MoSe chains with a 
lattice spacing clq = 0.85 nm, separated by Li atoms 
(Fig. n^a)). Atomic scale bundles of MoSe nanowires 
were produced from ^100 /LtM solutions of Li2MoeSe6 in 
anhydrous methanol. The solutions were then spun onto 
degenerately doped Si/SiG-2 substrates with lithograph- 
ically patterned electrodes in a Nitrogen atmosphere. 
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FIG. 1: (color online) (a) Structural model of a 7-chain MoSe 
nanowire along with the triangular Mo3Se3 unit cell, (b) and 
(c) AFM height images of MoSe nanowires between two Au 
electrodes. The wire heights are 7.2 nm and 12.0 nm respec- 
tively. Scale bar = 500 nm. 



Typically, 35 nm thick Au electrodes with 5 nm Cr ad- 
hesion layer separated by ~ 1 fim were used to contact 
randomly deposited nanowires. Figs, ^b) and (c) show 
atomic force microscope (AFM) images of typical de- 
vices. The two-probe resistance of such a device, which 
ranged from ~ 100 kil— 100 MCI at room temperature, 
was measured in a cryostat with a continuous flow of 
helium. A degenerately doped silicon substrate, under- 
neath a L ox — 300 nm thick silicon dioxide dielectric 
layer, served as a back gate to modulate the charge den- 
sity in the nanowires. Once transport measurements were 
complete, the wire diameter, D, was determined from an 
AFM height profile. 

Fig. |21 shows the conductance (G) normalized by its 
room temperature value as a function of temperature 
(T) for a representative subset of the samples studied 
[2fj . We applied a small bias voltage, (V <C fc^T/e), to 
stay in the linear response regime for this measurement. 
Notably, the mesoscopic scale samples (D < 20 nm, 
L <~ 1 /im) exhibit more than two orders of magnitude 
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FIG. 2: Relative conductance (G/G300K) versus T for six 
mesoscopic and two bulk MoSe wires labelled with the wire 
diameter. Curves are offset vertically for clarity. The solid 
lines are power-law fits to the data (G ~ T a ). The dashed 
line indicates the region where Coulomb blockade becomes 
important in two wires (A and T). Inset: Gate dependence 
of conductance (G) for the 5.3 nm wire (a) from 4K to 10 K 
at IK intervals. 

conductance decrease with decreasing temperature in the 
measured temperature range, unlike the samples in the 
bulk limit (D > 1 /im, L ~ 100 /im), which exhibit a 
bulk metallic behavior, as the conductance increases with 
decreasing temperature, i.e. energy. A power-law depen- 
dence, G ~ T a , is evident in these mesoscopic samples, 
where the exponent a can be readily extracted from the 
slope of the least-squares fit line in the double logarithm 
plot. For most of the samples, G(T) can be expressed 
by a single a within the experimentally accessible con- 
ductance range j2l| . However, for some wires with a 
relatively high conductance (> 1 (x$) at room tempera- 
ture (Fig. |21 ▲ and T), an abrupt change in the exponent 
at low temperatures was observed. In this low tempera- 
ture regime, the conductance varied with the gate voltage 
(Fig- HI inset) and the exponent depended on the applied 
gate voltage, V g , due to Coulomb charging effects. Above 
this Coulomb charging temperature, a general trend of 
decreasing a with increasing D is found for all meso- 
scopic samples studied. This trend will be discussed fur- 
ther later in the paper. 

We now consider several possible explanations for see- 
ing a decreasing conductance with decreasing tempera- 
ture. For example, one expects ln{G) cx — 1/T for bar- 
rier activated transport. For highly defective wires, one 
expects ln(G) oc —1/T 6 due to variable range hopping 
between localized states in the wire, where 5 can range 
from 1/4 for a 3D wire to 1/2 for a ID wire . However, 
neither of these models fit our data as, for all mesoscopic 
scale measured irrespective of wire diameter, the conduc- 



tance follows a power-law remarkably well. Another po- 
tential explanation is a non-metallic behavior associated 
with a Peierls transition, which opens up an energy gap 
at the Fermi level. However, from the conduction mea- 
surements on the bulk quasi- ID crystals (D > 1/im) and 
also from previous scanning tunneling microscopy work 
on similar nanowires we observe no evidence for a 
gap opening at temperatures down to 5 K, consistent 
with band calculations ^^|. Alternatively, a power law 
dependent tunneling conductance is predicted for tun- 
neling into a Luttinger Liquid, for ID Wigner Crystals 
[lif . for a highly disordered systems where the electron 
mean free path is comparable with the wire diameter 
We can rule out a ID Wigner Crystal model since 
the Coulomb interactions, which are screened by the back 
gate, are not long-ranged. We also eliminate the scenario 
for a strongly disordered system by estimating the elec- 
tron mean free path, l e in the nanowire. We can estimate 
l e indirectly from the effective wire length, L e tf, obtained 
from the dependence of the Coulomb charging energy on 
the gate voltage. Here, the charging energy E c — e 2 /C, 
where the wire capacitance C « 2ireL e ff /ln(4L ox / D). 
For the nanowire device shown in the inset of Fig. |2 the 
estimated charging energy E c is 5-10 meV from the con- 
ductance map in V and V g (not shown), from which we 
obtained L e ff ~ 0.3 — 0.6 /im. With this estimate, we 
rule out the possibility of having a highly disordered sys- 
tem. Moreover, the fact that the measured resistance is 
larger than ^100 kfl and is not directly correlated with 
the wire diameter indicates that transport is dominated 
by tunnelling. Thus a model concerning tunnelling into 
a relatively clean LL is a more likely description of the 
observed transport phenomena |22|. 

Further support for the LL-like transport in the MoSe 
nanowires can be found in the bias dependence of the 
conductance in the non-linear response regime. Accord- 
ing to the LL model in a tunneling regime the bias 
voltage dependent transport current, I(V), has a tran- 
sition between an Ohmic behavior, i.e. I cx V in the 
low bias regime (V -C fcsT/e), and a power law behav- 
ior with an exponent (3, i.e., / oc V@ +1 in the high bias 
regime (V ^> /csT/e). The inset of Fig. 3 shows typical 
I{V) data measured in a mesoscopic wire with the ap- 
plied bias voltage ranging over more than three orders 
of magnitude at different temperatures. A transition be- 
tween Ohmic and power-law behavior is observed as V 
increases. 

Interestingly, we also found that the exponent f3 de- 
pended strongly on D, as can be seen in Table [fl where 
we list a, (3 and D for 13 samples. In general, we find that 
a decreases monotonically as D increases. Based on the 
relation between a and (3, we can categorize our samples 
into two distinct groups: group (I) where a « 2/J; and 
group (II) where a « (3. In our experiments, the majority 
of samples (10 out of 13 in Tabled belongs to group (I), 
while only a few samples (3 out of 13 in Table HJ belong 
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TABLE I: Measured exponents a and /3 determined from the temperature and bias dependent conductance measurement, along 
with the wire diameter (D) as determined by AFM and the number of channels including spin (N) calculated from D. Wires 
indicated by asterisk (*) have a m f3 but for all other wires, a ~ 2/9. 



Wire 


Wl 


W2* 


W3 


W4 


W5 


W6 


W7 


W8* 


W9 


W10 


Wll 


W12* 


W13 


a 


6.6 


5.2 


4.3 


3.95 


2.33 


1.40 


2.34 


1.1 


1.55 


1.95 


1.2 


0.94 


0.61 





3.0 


4.9 


2.1 


1.9 


1.0 


0.72 


1.2 


1.0 


0.8 


1.09 


0.6 


0.90 


0.32 


D 


0.8 ± 0.5 


2.1 ± 0.3 


3.0 ± 0.3 


3.5 ± 0.2 


5.0 ± 0.5 


5.3 ± 1.0 


6.1 ± 0.5 


7.2 ± 0.5 


7.3 ± 1.4 


7.4 ± 1.5 


10.3 ± 0.4 


12.0 ± 0.7 


15.7 ± 1.1 


N 


2 


12 


22 


30 


62 


70 


94 


130 


134 


138 


268 


362 


620 



to group (II) 23]. 

For a clean LL without defects, these two exponents 
are expected to be identical, i.e., a = f3 Q, since they 
are characteristic of a single junction between FL and LL. 
The deviation from this model found in group (I) samples 
can be explained within the LL model with a few defects 
as described below. Strong defects in the nanowires break 
the conducting channels into a few LL dots connected in 
series between the electrodes. In this multiple LL dot 
scheme, the wires have two kinds of tunnel junctions; (i) 
junctions between the electrode and wire, constituting a 
FL to the end of LL junction; and (ii) junctions between 
two wire segments, constituting an LL to LL junction. 
In such a wire, the tunneling probability can be specified 
by two distinct exponents a_y,_LL and otpL-LL, where 
a ll—ll = 2a fl—ll holds |3J. If the linear response 
resistances of the FL-LL and LL-LL junctions at room 
temperature are of similar orders of magnitude, we ex- 
pect G tx T aLL - LL in the low temperature limit, since the 
LL-LL junctions become most resistive. However, in the 



high bias limit the FL-LL junctions become more resis- 
tive than the LL-LL junctions since (Xll-ll > o>fl-ll, 
and thus / oc v aFL ^ LL+1 . Therefore the exponents ob- 
tained from temperature and bias dependent data are 
expected to have a relation a = 2(3 for wires with a few 
defects that break them into multiple LL dots, as ob- 
served in our group (I) samples [3| • This argument does 
not hold however if there are no defects in the wire or 
in the extreme of strong defects that dominate transport 
within the experimentally accessible range of T and V. 
For such samples, a « (3, which is the relation we find in 
group (II) samples. 

The power-law behaviors in T and V allow us to scale 
I(V, T) into a single curve [3tll8|, Considering the above 
arguments, we can modify the scaling formula of a clean 
LL transport model to include the two exponents a and 
(3 as: 
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FIG. 3: (color) Inset: Wire W3 I-V data taken at differ- 
ent temperatures between 35 K and 200 K. All curves show 
a change from linear response to power-law dependence at a 
temperature dependent bias voltage. Main panel: I/T a+1 de- 
termined from I-V data plotted against jeV/ksT with Eq. 
fit to the data. The measured exponents are a = 4.3 and 
j3 — 2.1 (~ a/2) and the fitting parameter 7 is 0.25 ± 0.1. 



where Iq and 7 are constants independent of T and V. 
Physically, 7 represent the ratio between the voltage 
across a dominantly resistive junction at high bias to the 
applied bias voltage |3J. As shown in Fig. [31 the series of 
I(V) curves measured at different temperatures for the 
same sample collapse remarkably well onto a single curve 
described by Eq.^over the entire measured temperature 
range by plotting J /T a+1 against jeV/ksT with only 
one fitting parameter 7. For the data in Fig. 7 is 
0.25 ±0.1, implying that there are probably four barriers 
of approximately equal resistance over which the applied 
bias voltage is distributed. 

Finally, we now discuss the dependence of the expo- 
nents on the wire diameter (D) in order to elucidate the 
transition from a few channel ID transport to the 3D 
transport limit. For this purpose, we focus on the sam- 
ples with a — ctLL-LL (i-e. group (I) samples and the 
group (II) samples with 7 ss 1). In Fig. 21 we show the 
measured a plotted against D. Since N cx D 2 , where N 
is the number of channels in the wire including the spin 
degree of freedom |25|. the observed rapid decrease of a 
for large diameter nanowirc bundles indicates a cross over 
from ID behavior to 3D transport (a « 0). Employing 
the electron-electron interaction parameter for a single 
chain (N=2), g, the end tunneling exponent olll-ll for 
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FIG. 4: The exponent a, plotted against the wire diameter 
D. The solid line is a fit to Eq.E] 



an N channel LL wire can be expressed as 0, E| : 



olll-ll 



2 r 

N 



(1 + NU) 1/2 - 1 



(2) 



where U ~ 2/g 2 . We fit this equation to our data using 
g as a single fitting parameter. A good agreement with 
the experimental observation was obtained for g = 0.15. 
For a screened Coulomb interaction, a can be estimated 
by g ~ y/ \ / (e 2 ln(AL ox / aa) / irhv f k) where vf is the 
Fermi velocity of a single chain and k is the dielectric con- 
stant of silicon dioxide. From the fit in Fig. 21 we deduce 
that »f=3x 10 4 m/s. We note here that this value is 
smaller than the value obtained from recent band calcu- 
lation (4 x 10 5 m/s) jlif. indicating that a static screening 
picture considered in this model might be too simplistic. 
Further theoretical considerations including the effect of 
impurities and inter-chain hopping are needed to eluci- 
date strongly interacting electrons in these ID channels. 
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The mesoscopic wires were susceptible to structural de- 
formations locally at the wire/electrode junctions. Due 
to the invasiveness of the electrodes, multi-terminal mea- 
surements were not possible. 

Since G decreases rapidly as T decreases in the nanowire 
samples, most of our G(T) data were limited at low T by 
our current sensitivity (10~ 13 yl). 

We note that the environmental Coulomb blockade 
(ECB) model (see for example, Ingold and Nazarov, 
cond-mat/0508728 ) might be applied to our experiment. 
However, its validity is limited to the large N limit |8j, 
where the LL and ECB models do not offer any observ- 
able differences in our experimental setup and we thus 
consider only the LL model in this paper. 
Group (II) samples can be sub-divided further into a 
group with notable defects such as the one shown in 
Fig. 0[b) (7 ~1) and a group without strong defects 
(7 ml/2), where the parameter 7 is defined in Eq.0 
For temperatures below 300 K, the deviation of the 
measured a from oll-ll is less than 10% as long 
as R L ead-wire/ R-Defect is between ~ 0.1 - 10 where 
R-Lead-wire is the total resistance between the wire and 
the Au lead and RDefect is the total resistance of all de- 
fects. A similar tolerance range for Rhead-wire/ RDefect 
is also found for ~ oifl-ll- We have observed a few 
devices outside of this tolerance limit, which shows an 
intermediate bias ranges where LL-LL tunnelling is dom- 
inant at low temperature. 

We use N = 2fD 2 /al where /, the filling factor, is 0.91 
assuming a hexagonal close packing of the chains. 



